From chaos to disorder: Statistics of the eigenfunctions of microwave cavities 
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We study the statistics of the experimental eigenfunctions of chaotic and disordered microwave 
billiards in terms of the moments of their spatial distributions, such as the inverse participation 
ratio (IPR) and density-density auto-correlation. A path from chaos to disorder is described in 
terms of increasing IPR. In the chaotic, ballistic limit, the data correspond well with universal 
results from random matrix theory. Deviations from universal distributions are observed due to 
disorder induced localization, and for the weakly disordered case the data are well-described by 
including finite conductance and mean free path contributions in the framework of nonlinear sigma 
models of supersymetry. 
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I. INTRODUCTION 

This paper briefly summarizes results and insights 
gained from expeirmental studies of eigenfunctions of 
chaotic and disordered 2D billiards. The experiments 
which utilize microwave cavities, enable us to tune the 
degree of localization by varying two parameters, the fre- 
quency / and the mean free path I. Thus we are able to 
access the ballistic limit I A (the wavelength) as well 
as the strongly localized limit I -< A. In our earlier work 
jjj we have experimentally shown that the spatial inten- 
sity distribution of an eigcnfunction of a chaotic billiard 
quantitatively agrees with the Random Matrix Theory 
(RMT). Recently, we have obtained further results on 
the statistics of the eigenfunctions || of chaotic and dis- 
ordered systems, which are reviewed in this paper. The 
results of our experiments on disordered microwave cav- 
ities are quantitatively described by calculations based 
upon the supersymmetry sigma models carried out by 
Efetov ||, Prigodin and Altshuler 0, Mirlin g and oth- 
ers. In this paper we explore the eigenfunction statistics 
of chaotic and disordered systems using inverse partici- 
pation ratio as a measure of the disorder strength. 



II. THEORETICAL BACKGROUND: RMT AND 
NONLINEAR SIGMA MODELS 

It was shown by Berry that the spatial amplitude dis- 
tribution of the wave-function of an electron in a chaotic 
cavity can be generated from the random superposition 
of plane waves. This leads to the universal spatial distri- 
bution of the eigenfunction P(^>) which is Gaussian, and 
which implies that the spatial intensity distribution of 
the eigen function Pd^l 2 ) is Porter-Thomas ||. Spatial 
correlations are also obtained, such as the density-density 
autocorrelation (\^(r)\ 2 \^(r')\ 2 \ = l + 2J 2 (k\r -r'\) for 
& 2 — D chaotic billiard. 



The moments of the spatial integral of density of an 
eigenfunction /„ = J \^l(f)\ 2n d 3 r are important mea- 
sures of the localization, particularly the second moment, 
the Inverse Participation Ratio 1% (IPR), and its statis- 
tics P/ 2 (/2), which describes important properties of the 
chaotic and disordered systems. The eigenfunctions of 
a classically chaotic billiard are delocalized wave func- 
tions in the limit of an infinite system. The moments I n 
have fixed values with no fluctuations, with I2 = 3.0 in 
2-D, i.e. Pi 2 {h) = 6(I 2 - 3.0) @§. The universalities 
have been studied for some time within the theoretical 
framework of the RMT. 

From the perpective of transport in disordered sys- 
tems, the chaotic system can be considered as the bal- 
listic limit corresponding to infinite mean free path and 
infinite conductance. Due to the presence of weak dis- 
order, the universal properties break down, and system 
properties deviate from the universal values. A perturba- 
tive treatment has been achieved using nonlinear sigma 
models of supersymmetry, originally motivated by the 
problem of electrons in a disordered metal, which can be 
termed quantum diffusion (?]]. The supersymmetry ap- 
proach, which treats an electron in a potential Vi mp (r) 
characterized by a Gaussian distribution of random disor- 
der due to impurities, enables the calculation of spectral 
correlations, such as level spacing statistics, form factor, 
etc. and eigenfunction correlations including amplitude 
P(^>) and density P(|<3>| 2 ) distributions, and spatial auto- 
correlations (||^'(r)| 2 |^'(r )| 2 Y In the spatially homoge- 
neous 0-mode limit of the theory, which corresponds to 
the ballistic limit kl » 1, the theory reduces to that of 
the chaotic limit and yields the results of RMT, such as 
the Porter-Thomas density distribution, characteristic of 
systems with a Gaussian distribution of amplitudes. 

With increasing disorder, i.e. with a finite mean free 
path and a finite conductance, the system statistics devi- 
ate from the RMT, and now the system can be described 
by ID sigma models as a perturbative correction with a 
finite conductance g and a finite disorder controlled pa- 
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rameter kl, where k is the wave vector and I is the finite 
mean free path. The leading term, the 1-mode limit, 
yields the leading corrections due to incipient localiza- 
tion, with higher orders (in principle) leading to increas- 
ing localization. The intensity-intensity auto-correlation 

| x P(r)| 2 |^(r')| 2 ^) also can be described by sigma model 

calculations. For a disordered system, correlation is large 
at r — r' = and dies out as with the increasing distance. 

III. EXPERIMENTAL DETAILS 

Recent developments in experimental techniques have 
made it possible to measure the eigenvalue and eigen- 
function properties of thin cylindrical microwave cavities 
- these studies have enabled quantitative studies of issues 
in Quantum Chaos and Localized media on laboratory 
length scales with very high precision experiments. 

In 2D, Schrodinger and Maxwell equation map onto 
each other in the sense that the Hclmholtz equation for 
the z-component of the electric field becomes (V 2 + 
k 2 )ip — 0, where tp — E z . 2D cavities were made be- 
tween two parallel copper plates, and the shape and size 
of a cavity was made according to the need. The height 
of the cavity is d = 6 mm, and the cavity is effectively 
2D for the range of microwave wave lengths used in the 
experiments. The electric field distribution for the Trans- 
verse Magnetic (TM) modes of such a cavity obeys the 
Schrodinger wave equation in two dimensions |uj . 

For chaotic billiard experiments, the shape of the cav- 
ity was made like a Sinai Stadium or like Sinai Billiards. 
For disordered closed billiard experiments, the cavity has 
closed rectangular ( 44cm x 21.8cm ) boundaries, and 
the disordered scattering centers were imposed by placing 
randomly circular 2D plates of 1cm radius. The scatter- 
ing mean free path was changed by changing the number 
of the scattering centers. Several realizations of a scatter- 
ing system were done by changing the scattering positions 
with a random number generator and keeping the total 
number of scatterers as same. A variety of experiments 
on closed and open quantum chaotic systems have been 
carried out using such microwave cavities (s| |l0| . A par- 
ticularly powerful aspect of the experiments is the ability 
to directly measure wavefunctions using cavity perturba- 
tion techniques 

Experiments on the eigenvalue statistics have enabled 
stringent tests of the eigenvalue properties of chaotic sys- 
tems Q . Experiments on disordered billiards show a very 
good realization of 2D disordered systems , and which 
motivated further theoretical works in this field pj],^,p[. 
Most recently, we have carried out a systematic study 
of the statistics of the experimental eigenfunctions for 
chaotic and disordered systems gj. These papers show 
tests of several analytical results of RMT and the nonlin- 
ear model of supersymmetry which was developed mainly 



for disordered granular media and mesoscopic systems. 
At present, the eigenfunctions of mesoscopic systems are 
difficult to access for controlled experiments due to their 
smaller size, and consequently microwave cavity experi- 
ments play an unique and important role for wave func- 
tion statistics. 

IV. EXPERIMENTAL EIGENFUNCTIONS AND 
THEIR STATISTICS 

The utility of analyzing eigenfunctions in terms of their 
IPR is discussed here for a few representative eigen- 
functions. For convenience, we will use the notation 
h=J |*(r)| 4 dw = Ju 2 dv, u = |*(r)| 2 , w = (I 2 - 3)/6, 
and I\= J ^(r)) 2 ^ = 1, where the integration is done 
over the 2D volume. More than 250 eigenfunctions are 
analyzed for chaotic and disordered systems, and each 
eigenfunction has 3200 spatial data points. 

The representative eigenfunctions of the disordered 
billiard with N = 36 scatterers have been shown in 
Figs. 1(a), (b), (c) and (d). Fig. 1(a) shows the eigenfunc- 
tion of a disordered billiard with a low eigenfrequency 
/ = 3.84GHz has localized eigenstate IPR with value 
I 2 = 11.22. With increased frequency (energy), the same 
disordered medium shows more delocalized eigenstatcs 
with decreasing IPR values : Fig. 1(b) I2 = 9.64 at 
4.49GHz, Fig.l(c) I 2 = 7.52 at 6.37GHz, and Fig.l(d) 
IPR value I 2 = 4.36 at 7MGHz. This can be regarded 
as a tuning the degree of localization by changing the 
frequency (energy). In fact, we have shown that this 
tuning path follows a power law decay with an exponent 
i , consistent with theoretical calculations Jl2| . 

Fig. 1(e) and (f) show representative eigenfunctions of 
a chaotic billiard, which are "delocalized", IPR I 2 = 
3.13 at 6.05GHz (Fig.l(e)), and I 2 = 3.02 at 7.54GHz 
(Fig. 1(f)). These eigenfunctions have their IPR values 
close to the universal value (I 2 ) = 3 in 2D. 

A comparison of experimental results of chaotic and 
localized eigenfunction statistics are shown in Figure 2 
for (i) spatial density distribution of eigenfunctions, (ii) 
statistics of eigenfunctions interms of their IPR values, 
and (iii) intensity auto-correlation of the eigenfunctions. 

A. Spatial density distribution of eigenfunctions 

Figure 2(a) shows that the spatial intensity distribu- 
tion Pd^l 2 ) of eigenfunctions of a chaotic system obeys 
the Porter-Thomas distribution. The P-T distribution 
is obtained from RMT, as well as random superposition 
of plane waves and the 0-D sigma model. The deviation 
from the PT distribution due to finite localization of a 
disordered medium is shown in Figure 2(b). This can 
be modeled by introducing finite mean free path in the 
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problem and a quantitative description of the data has 
been achieved 

B. Statistics of IPR of eigenfunctions 

The distribution P/ 2 (/ 2 ) of the IPR ^for the chaotic 
billiards is shown in the Fig. 2(c). The distribution is 
nearly symmetric around (I2) = 3. RMT calculations 
indicate that there should be no fluctuations, but this 
only applies for an infinite system. The finite width of 
the distribution is due to the boundary scattering, and 
can be quantitatively described by large but finite 2kl 
and g values [fi3| , A quantitative description is given 
later. 

In contrast the IPR distribution (Fig. 2(d)) for a dis- 
ordered billiard with mean free path I — 5.1cm (N = 36 
scatterers) is strongly asymmetric and has a mean value 
much larger than the universal value of 3.0. These data 
can be modeled with a finite conductance g . A dimen- 
sionless conductance can be defined as g — ln(R/l)/{w), 
where R is the system size and / is the mean free path 
and (..} is the realization average for a fixed "disordered 
strength" 2kl. 

Recent theoretical calculations by Prigodin and 
Altschuler have shown that when h > (h) , the IPR 
distribution follows an exponential decay law 

Pj a (/ 2 ) = c72i/f eM—ah) (1) 

y I2 6 
while P{I 2 ) for h < (h) is § 

PiAh) = Cx|exp[-|(7 2 - (J 2 » - I e -*a»-<&»], (2) 

where C\ and C 2 are the normalization constants. For 

the sample with finite mean free path I = 5.1cm, and 
assuming random phase approximation within a small 
window where 2kl is fixed and independent of IPR dis- 
tribution [fL8| , experimental data matches well with Eq.l 
with a conductance value g = 1.0 as shown in Fig. 2(d). 
The distribution is also asymmetric beyond the universal 
average value (I2) = 3, for the chaotic billiards. 

For the chaotic billiard, the data matches well with 
Eq.2, as shown in Fig. 2(c) for the finite conductance 
g = 7.8 and 2kl — 37, quite large values that we can 
consider the system as infinite, i.e., the length scale same 
as the system size. Values of g are consistent with the 
parameters of the experimental microwave cavity. 

C. Intensity-intensity auto-correlation of the 
eigenfunctions 

In earlier work we have shown that the intensity- 
intensity auto-correlations for the chaotic cavity has a 



universal behavior associated with Friedel oscillations 
Jl9| . We have extended the density auto-correlation cal- 
culation from chaotic to localized eigenfunctions, and 
analyzed the results analytically and numerically. The 
intensity-intensity auto-correlation for an arbitrary dis- 
ordered strength 2kl can be calculated as follows. Let us 
define K(r) = | ImG(r')\ 2 /(ttv) 2 , where G(\r-r\) = 
(r\(E — H)~ 1 \r 1 ') is the Green function of the disordered 
system Hamiltonian, then one can show that |2(]] 
K(r) = ^lZoT^Mkr(l + ^y)]dy\ 2 ._ 
where Jo is the zeroth order Bessel function. 
For chaotic billiards the auto-correlation follows an an- 
alytical form 

{^(r^ 2 (r')) = l + (I 2 -l)J 2 (r-r') (3) 

This also corresponds to that for a fiat disordered po- 
tential with Gaussian fluctuations, when 2kl » 1. We 
extend the previous calculations for the finite localization 
case. Repeating the earlier calculations of |l^] with a fi- 
nite kl value, it can be easily shown that, the expression 
for the auto correlation can be approximately expressed 
in a region r — r' ^ I by a decay length scale of scattering 
mean free path length I as follows. 

(tf 2 (r)*V)) ~ 1 + {h - l)K(k\r - r'\) (4) 

~l + (/ 2 -l)J 2 (fc|r-r'|)e-^^ (5) 

With kl large, the above expression, Eq.4 and 5, will 
converge to the expression of a chaotic system, Eq.3. 

We have plotted the value of (^> 2 {r)^ 2 {r')) for differ- 
ent types of chaotic billiards in Fig. 2(e) and it agrees 
well with Eq.3. 

For the disordered billiards, K was solved numerically, 
and also approximate analytical results have been plot- 
ted in Fig. 2(f), and were compared with the results of 
(^' 2 (r)^I' 2 (r')) calculated from experimental data. The 
agreement is quantitatively excellent in the described re- 
gion. For Sinai billiards, 2kl = 37 is a very large num- 
ber and described well for a Gaussian distribution of the 
wave functions (Fig. 2(f)). For the disordered billiards, 
the matchings are shown for the 2kl = 7 and 3 (Fig. 
2(f)). The correlation is large at r — r' = 0, i.e. equals 
(\I/ 4 (r)) = I 2 , and the correlation decays to zero with in- 
creasing spatial distance, as the auto-correlation is neg- 
ligible at a larger distance beyond the scale of the local- 
ization length. 

V. CONCLUSIONS 

In this paper we have reviewed statistical properties 
of wave functions of chaotic and disordered systems us- 
ing the second moment of the spatially integrated in- 
tensity of the eigenfunctions. We have shown that 
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the statistics agrees well with random matrix theory 
for chaotic systems, and with the leading perturbative 
correction of the nonlinear sigma model of super sym- 
metry for incipient disordered media. The perspective 
that emerges from the microwave experiments is summa- 
rized in Figure 2 in terms of spatial intensity distribu- 
tion (Fig.2(a,b)), statistics of Inverse Participation Ra- 
tio (I PR) of eigenfunctions (Fig.2(c,d)) , and intensity- 
intensity auto-correlation function (Fig.2(e,f)). The uni- 
versal (2D) results of Quantum Chaos, characterized by 
Gaussian fluctuations of are achieved in chaotic bil- 
liards, albeit with small corrections due to correlations 
induced by boundary scattering. The importance of cor- 
relations due to coherent interference between scattering 
events are amplified in the model disordered billiards. By 
tuning kl we can observe a path from extended (1% ~ 4) 
states to strongly localized (I2 ~ 20) states, as can be 
seen in a I2 vs. f 2 plot (Fig. 2 in Ref. Q). This is sim- 
ilar to tuning from a metal to insulator by increasing 
energy! This is also accompanied by deviations from 
the Porter-Thomas (PT) distributions, and strong level- 
to-level fluctuations, leading to a strongly asymmetric 
distribution Pj 2 (/2), and strong spatial decay of corre- 
lations, all indicating non-Gaussian amplitude distribu- 
tions. For not too strongly localized states, we have 
shown agreement with the supersymmetry theories, pa- 
rameterized by finite conductivity g (Figures: 2(b), 2(d) 
and 2(f) ). The conductivity g depends monotonically 
on kl, with the ballistic chaotic limit corresponding to 
g, kl — ► 00 . A systematic study of wave function statis- 
tics is under way to explore more strongly localized states 
for which however there is no available theory currently. 

The eigenvalue and eigenfunction statistics of a dis- 
ordered medium for finite g have been calculated Q by 
solving a Perron-Frobenius operator of a classical diffu- 
sion equation, whose (real) eigenvalues are related to the 
poles of a classical zeta function. Recently, a direct ex- 
perimental observation of the complex poles of the Liou- 
ville operator for an open chaotic system, the so-called 
classical Ruelle-Pollicot resonances, have been observed 
by us in microwave experiments on open n-disk systems 
|i"4]-|i"7f . This latter work thus provides a physical reality 
to these eigenvalues in open systems which have played 
a significant role in the theory of the disordered systems, 
and provides a close connection between classical and 
quantum diffusion. 

While the nonlinear sigma model provides a useful, 
even quantitative, framework to discuss arising from co- 
herent interference leading to localization, its perturba- 
tive approach limits it to the leading corrections to the 
ballistic or RMT limit. The experiments are not lim- 
ited to the regime of large g i.e. kl ^> 1, but can easily 
reach the strong localization limit kl ~< 1. This regime re- 
quires higher order calculations beyond the 1-mode that 
are quite formidable. Other approaches to this problem 



deserve attention. One possibility that intrigues us is the 
information theoretic approach based upon maximum en- 



tropy employed by Kumar and collaborators |21|. While 
this has thus far been applied only to quasi ID, exten- 
sions to higher dimensions would be worth pursuing. 

It is a pleasure to dedicate this paper to Professor 
Narendra Kumar. PP had the honor of carrying out 
his Ph.D. thesis research under Professor Kumar's guid- 
ance. SS has had the pleasure of numerous discussions 
on physics that have always been enthralling and memo- 
rable. This work was supported by NSF-PHY-0098801. 
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FIG. 1. (a,b,c,d) Representative eigenfunctions of a disor- 
dered billiard with number of scatterers is 36 : (a) IPR value 
h = 11.22 at / = 3.84GHz, (b) J 2 = 9.64 at 4.49GHz, (c) 
h = 7.52 at 6.37GHz and (d) h = 4.36 at 7.20GHz. The 
plots show that the eigenfunctions are more localized ( larger 
IPR value) at a lower frequency and more delocalized (smaller 
IPR value) at a higher frequency. 

Fig.l(e,f) Representative eigenfunctions of a chaotic billiard: 
(e) h = 3.13 at 6.05GHz , and (f) h = 3.02 at 7.54GHz. 

FIG. 2. (a,b) . The spatial intensity distribution of exper- 
imental eigenfunctions: 2(a) Porter- Thomas distribution for 
chaotic cavity, and 2(b) non Porter-Thomas distribution for 
disordered cavity. 

Fig. 2(c,d). 2(c) The IPR Pi 2 (I 2 ) distribution of the 
chaotic Sinai-stadium billiard is nearly symmetric around 
the universal mean value 3.0 . 2(d) The IPR distribution 
Pr 2 (/2) of the disordered billiards is strongly asymmetric and 
non-Gaussian. The lines represent calculations based on the 
nonlinear sigma-models. 

Fig 2(e,f). Intensity auto-correlation (^ 2 (r)^f 2 (r')^ of eigen- 
functions of a 2(e) Sinai-stadium billiards, and 2(f) disordered 
billiards with different fixed disordered strengths 2kl, exper- 
iment ( dotted lines), Eq.(4) (dashed lines), and numerical 
solution of Eq.(4) with 2kl =13 and 7 (solid lines). 
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